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ABSTRACT 

Aims. We extend the models presented in Mordasini et al. (2009) to the formation of planets orbiting stars of different masses. We 
discuss the properties of the resulting synthetic planet population in terms of mass, orbit, and metallicity distributions. 
Methods. The population synthesis calculations we use are based on the planet formation model developed in Alibert et al. (2005a), 
which self-consistently takes into account planetary growth and migration in an evolving proto-planetary disk. Using this model, we 
generate synthetic populations of planets by following their growth in a large number of proto-planetary disks, whose properties (mass 
and lifetime) are selected in a Monte Carlo fashion using probability distributions derived from observations. 

Results. We show that the scaling of the proto-planetary disk mass with the mass of the central star has a direct and large influence 
on the properties of the resulting planet population. In particular, the observed paucity of large mass planets orbiting 0.5 Mq stars can 
directly be explained as resulting from a only slightly steeper than linear scaling. The observed lack of short period planet orbiting 2.0 
Mq stars can also be attributed to this scaling but only if associated with a decrease of the mean disk lifetime for stars more massive 
than 1.5 Mq. Finally, we show that the distribution of minimum mass and semi-major axis of our synthetic planets are statistically 
comparable with observations. 

Key words. Stars: planetary systems - Stars: planetary systems: formation - Stars: planetary systems: protoplanetary disks - Planets 
and satellites: formation - Methods: numerical 



1. Introduction 

Since the pioneering discovery of the first extrasolar planet or- 
biting a main sequence star by Mayor and Queloz (1995), more 
than 400 extrasolar planets have been discovered, the majority 
of them by means of the radial velocity technique. After having 
essentially focussed on solar type stars, more and more surveys 
are targeting stars of different nature, in particular of different 
mass. From a theoretical point of view, these surveys are impor- 
tant since a global understanding of the processes at work during 
planet formation relies on the comparison of observations with 
model results in environment as varied as possible. 

In the core-accretion scenario of planet formation, gas gi- 
ant planets are believed to form in two steps. In a first step, 
the planet's core is formed by accumulation of planetesimals, 
in a way similar to the one generally accepted for the forma- 
tion of terrestrial planets. When a critical mass is reached, the 
core starts to capture ever increasing quantities of nebular gas, 
leading to a phase of runaway gas accretion at the end of which 
the planet reaches its final mass. For a long time, this forma- 
tion model suffered from a nagging timescale problem. Unless 
the disk was quite massive, the formation time, as computed by 
the model, was found to significantly exceed observationally in- 
ferred proto-planetary disk lifetimes. In recent years, improv- 
ing the physics handled by the models and/or considering ad- 
ditional physical processes has led to a significant speed-up of 
the formation process thereby removing this timescale problem. 
Indeed, theoretical models based on new opacity determinations 



of Podolak (2003) allowed formation of Jupiter-like planets in 
a few Myr (see e.g. Hubickyj et al. 2005). Moreover, extended 
core-accretion models taking into account migration of form- 
ing planets, as well as proto-planetary disk structure and evo- 
lution led to equally short formation time (see Alibert et al., 
2004, 2005a). Finally, these latter models have been shown to 
be compatible with our present day knowledge of Jupiters and 
Saturn's internal structure and atmospheric composition (Alibert 
et al. 2005b). 

Many different studies have already considered the forma- 
tion of planets orbiting stars of different masses. Laughlin et al. 
(2004) have shown for the first time that the formation of gas 
giant planets orbiting low mass stars is inhibited. This is mainly 
due to the reduced Keplerian frequency in this environmenj^ and 
to their assumption of a strong scaling between disk and primary 
mass. Ida and Lin (2005) used their planet formation model to 
study the influence of the central star mass on the characteristics 
of the planets. To relate the mass of the proto-planetary disk to 
the one of the primary star, they used a scaling relation of the 
type: Mdisk ^"tai-^ where ao was allowed to vary between 
and 2. Their study shows that close-in Neptune-mass ice-giant 
planets should be common around low mass stars, while Jupiter- 



' The Keplerian frequency Q = y "^'^'f^ , where Mstar is the mass of 
the primary, G the gravitational constant, and r the distance from the 
central star, gives the typical timescale of coUisional processes involved 
in the formation of planets. 
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mass gas-giant planets should be quite rare. On the opposite, the 
high-mass peak of the planetary mass function grows in impor- 
tance around high mass stars. 

Finally, Kornet et al. (2006) have calculated the influence 
of the central star mass on the evolution of the protoplanetary 
disk, by focusing in particular on the redistribution of solids in 
the disk. Using a simple model to relate the resulting structure 
of the proto-planetary disk and the probability of presence of 
a giant planet, they have concluded that, globally, the probabil- 
ity to harbour a Jupiter mass planet increases with a decreasing 
primary mass. Interestingly enough, their conclusions are very 
different from the ones by Laughlin et al. (2004) and Ida & Lin 
(2005). Indeed, in these studies the surface density of planetesi- 
mals increases with the mass of the primary, whereas in Kornet 
et al. (2006) it decreases due to the fact that the redistribution of 
soUds is more efficient in low mass star environments. Present 
day observations suggest that the probability to harbour an ob- 
servable giant planet increases with the mass of the central star, 
which shows that the redistribution of solids in protoplanetary 
disks is not the dominant process in the early stages of planet 
formation. 

In a recent paper, Mordasini et al. (2009a, paper I) have pre- 
sented planet population synthesis calculations based on the ex- 
tended core-accretion formation model of Alibert et al. (2005a) 
which takes into account migration and disk evolution. This 
model has been shown to reproduce some of the bulk proper- 
ties of the giant planets in our own Solar System (Alibert et al. 
2005b), and of the three Neptune mass planet system discov- 
ered around HD69830 (Lovis et al. 2006, Alibert et al. 2006). 
Population synthesis calculations rely on the calculation of the 
formation of thousands of planets, each one assuming a differ- 
ent set of initial conditions. In particular, the mass of the proto- 
planetary disk, its lifetime, its metallicity (more exactly the 
planetesimals-to-gas ratio in the disk), and the starting location 
of the planet's initial embryo have to be specified. The distribu- 
tion functions for these initial conditions are taken, as far as pos- 
sible, directly from astronomical observations. Furthermore, by 
taking into account the detection bias introduced by radial veloc- 
ity surveys (Naef, 2004), it was possible to compare in a statis- 
tical way the observed population of giant planets with the syn- 
thetic ones. Using Kolmogorov-Smirnov (KS) tests, Mordasini 
et al. (2009b, paper 11) showed that this approach yields a planet 
population whose mass and semi-major axis distributions match 
the properties of the actually observed planets to a confidence 
level of close to 90 %. 

In this paper, we extend the calculations presented in papers 
I and II to the formation of planets orbiting stars of different 
masses. In order to isolate the effects of varying the mass of the 
star, we use the same formation model varying only the central 
star mass. Note that, since we also consider models for which 
the properties of the proto-stellar disk vary with the mass of the 
primary, this also leads for these cases in changed disk charac- 
teristics. 

Since we not only use the exact same procedure as explained 
in detail in paper I and II but also all the same value for all nu- 
merical parameters, we wiU not discuss our approach in any de- 
tail here but refer the reader to these papers. The present paper 
is organized as follows: in Sect. 2, we discuss the influence of 
the mass of the primary on the formation process itself as well 
as on the corresponding initial conditions. In Sect. 3, the out- 
come of model in which the mass of the central star is varied 
from 0.5 Mq to 2.0 Mq are presented. In particular, we discuss 
the changes in the properties of the resulting planets (mass and 
semi-major axis, composition) as well as the influence of metal- 



licity in the formation process. Finally, Sect. 4 will be devoted to 
the discussion and conclusions. 



2. Formation model and initial conditions 

2.1. Influence of the central star mass 

The planet formation model we use is a simplified version of the 
one presented in Alibert et al. (2005a). As we have already de- 
scribed this model extensively in paper I, we will not repeat it 
here. Our approach consists in following the growth of a single 
planetary embryo initially located at a given distance of the cen- 
tral star in a disk consisting of gas and soUds. Note that the opac- 
ity is not reduced compared to interstellar values (e.g. Pollack et 
al. 1996, Hubickyj et al. 2005). In our calculation, the effect of 
such a reduction is indeed rather small. The reason for that is 
that thanks to migration (even with the strongly reduced type I 
migration rate), cores never get completely cut from a supply of 
fresh planetesimal to accrete, without the need of having the en- 
velope growing (where the opacity matters) in order to expand 
the solid feeding zone, as it is the case for a strict in situ forma- 
tion. These effects are discussed in Mordasini et al. (2010). 

The actual mass of the central star enters as a parameter in 
many of the physical processes involved during planet forma- 
tion. ExpUcitly, it enters in: 

(M 1 \^/-^ 
3^ j , where Mpianet and 

Mstar are the planet and star masses, and Opianet is the planet's 
semi-major axis). This radius is a measure of the size of the 
planet's feeding zone, and (indirectly) of its envelope size 
and hence mass. Since it inversely scales with stellar mass, 
planets forming in orbit of a more massive stars can accrete 
planetesimals originating from a smaller region of the proto- 
planetary disk. 

2. The type I migration rate. Type I migration rate, which is 
relevant for low mass planets, is reduced for more massive 
central stars (see Tanaka et al. 2002). 

3. The disk structure. On one hand, increasing the mass of the 
central star results in a stronger vertical component of the 
gravity in the disk. On the other hand, because the viscous 
dissipation rate is a function of Keplerian frequency, increas- 
ing the stellar mass leads to hotter disks. Numerical calcula- 
tions show that the first effect dominates over the second and 
that disks orbiting more massive stars, everything else being 
equal, are generally thinner. 

4. The gravitational radius Rg (see Veras and Armitage, 2004), 
which sets the region in the disk where gas can be photo- 
evaporated. This parameter depends linearly upon the central 
star mass (see Adams et al. 2004). 

5. The Keplerian frequency which, among other things, gov- 
erns, the accretion rate of soUds. 

In addition, there are stellar mass dependencies in the com- 
putation of the initial conditions. Initially, all our disks are 
characterized by a power law surface density of index -3/2. 
However, the temperature and pressure structures of these disks 
are functions of the mass of the central star Hence, the stellar 
mass determines the location of the iceline where the surface 
density of solids jumps by a factor of about four. All other pa- 
rameters being equal, the iceline moves outwards with increas- 
ing stellar mass by 1 to 2 AU (depending on the disk's mass) for 
a stellar mass varying from M^tar - 0.5Mq to Mstar - 2.0Mq. An 
analytical fit using our nominal value for the viscosity parameter 
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Fig. 1. Disk aspect ratio (thickness over distance to the star) as 
a function of distance to the central star. Each line represents a 
different epoch, from 0.1 Myr *uppermost line) to 5.2 Myr, with 
a 0. 1 Myr step. The initial mass of the disk is 0.016 Mq, and the 
lifetime is 5.3 Myr. The mass of the central star is 1.0 Mq. 

Q- = 7x 10"^ (Shakura & Sunyaev, 1973), provides the following 
approximation for the initial position of the iceline 

/ V \0-44 / ,^ \0.1 

^ice _ / ^i.lhV \ / ^"star | 

AU ~ \10g/cm2/ \ Mo / 

where Ss.iau is the initial gas surface density at 5.2 AU|^ 

We recall here that we only compute the disk structure in the 
inner 30 AU of the protoplanetary disk. This choice is governed 
by the fact that all planets, in these calculations, form in the cen- 
tral parts of the disk. In order to solve the diffusion equation, 
we have to assume a boundary condition at 30 AU. In the cal- 
culations presented here, we assume , as in Alibert et al. (2005), 
that the mass flux at 30 AU is null. In this sense, we overesti- 
mate the gaseous mass in the inner 30 AU, since part of the gas 
in reality goes beyond 30 AU for angular momentum conserva- 
tion. However, the viscous timescale at 30 AU is long enough 
so that the precise boundary condition there has a low influence 
on the planet forming region (below 10 AU or so). In addition, 
note that the lifetime of disk is, by construction, similar to the 
observed one, at least for solar type stars. Indeed, the photo- 
evaporation rate is adjusted so that our disk population has a 
disk lifetime probability consistent with observations. Note that 
changes in the disk lifetime only result from variations in the 
assumed photoevaporation rate and not from variations in the 
assumed viscosity parameter. 

Finally, as in Alibert et al. (2005), the structure of the disk is 
calculated without taking into account the irradiation of the cen- 
tral star, which itself would depend on the mass of the primary. 
Since higher mass stars tend to have higher luminosities, this 
would also translates into an iceline located at larger distances. 

In order to ease the comparison with other disk models. Fig. 
[T] presents the evolution of the disk aspect ratio (which is im- 
portant for the switch between type I and type II migration) for a 
"typical" disk, namely a disk with a mass of 0.016 M©, a lifetime 
of 5.3 Myr, around a 1.0 Mq star. 

Finally, the time required to grow the seed embryo used in 
our calculations depends on the Keplerian frequency as well as 

^ in our model, the location of the iceline is set by the initial condi- 
tions and does not vary with the evolution of the disk, see paper I 



Table 1. Parameters for the simulations used in the paper (except 
if expUcitly mentioned). 



Feature Values 

Type I migration reduction factor fi 0.001 

Viscous parameter a 7 x 10"' 

Initial exponent gas disk Z(a, f = 0) -3/2 

Rockline included no 

Iceline included yes 

Outer radius of the computational disk 30 

Inner radius of the computational disk 0. 1 

/d/g,o 0.04 



on the actual distribution of solids. This time determines the time 
offset between disk evolution and the start of embryo growth 
and hence the time available to accrete solids and gas before the 
disk disappears (see paper I). For a fixed surface density, numer- 
ical calculations show that the embryo formation timescale at a 
fixed distance decreases within increasing stellar mass. This ef- 
fect is further amplified by the fact that disks orbiting high mass 
stars are more massive than their counterparts orbiting low mass 
stars. We note that depending upon the location of the iceline, 
disks orbiting smaller mass stars may have a locally higher sur- 
face density than disks orbiting more massive stars, but this only 
concerns a small number of disks and is limited to small regions 
of these disks. 

2.2. Parameters and initial conditions 

In this paper, we aim at studying in a global framework the for- 
mation process of planets orbiting stars of different mass. As a 
consequence, we will use for the parameters of the model the 
same values as in paper I. Depending upon the scaling of the 
disk with the mass of the primary, the disk mass and lifetime 
will represent an exception to this rule. In addition, since part of 
our work focusses on the frequency of close-in planets predicted 
to orbit stars of various masses, we now follow the planetary mi- 
gration and growth until the planet crosses the inner boundary 
of the protoplanetary disk located at 0.1 AU. This means that 
the inner edge of the planet feeding zone, which was located at 
flpianet -4/?H in paper I and II, is now located at the greater of this 
former value, and 0.1 AU. This modification has no effect on the 
mass growth of planets. 

The two key parameters characterizing the formation model 
are the a parameter entering in the determination of the viscosity 
and the type I migration rate parameter //. As shown in paper II, 
population synthesis calculations allow to reproduce the semi- 
major axis and mass distributions for single giant planets orbit- 
ing G stars using a - 7x 10"^, and // - 0.001 . Other parameters 
of the model are taken from the nominal model in paper I, and 
are summarized in table [T] where /d/g.o is the assumed plan- 
etesimals to gas ratio in the planet forming region for a solar 
metallicity system. We stress that, for simplicity, we do not take 
into account the dependance of star's metallicity with respect to 
mass. Even if M dwarfs are known to be slightly metal poor, 
this effect is not included in our calculations. We justify this by 
noticing the much larger uncertainty existing on the relationship 
between disk mass and primary mass. 

The determination of the relation between the disk mass and 
the central star mass has been the subject of numerous obser- 
vational studies. In our model, we assume that the two quanti- 
ties are linked by a relation of the type Mjisk M"^^^. Based 
on Ha line profiles, MuzeroUe et al. (2003), and Natta et al. 
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(2004) have inferred that the disk accretion rate roughly scales 
with the square of the primary mass. Using our disk model we 
have determined the dependance of this accretion rate at the in- 
ner edge of our computational disk on the value of au- Our cal- 
culations show that ffD ~ 1.2 allows to roughly reproduce the 
data from Muzerolle et al. (2003) and Natta et al. (2004). We set 
this value as our nominal value but also consider the cases for 
which ao - and 2 in an effort to provide bracketing results. 
Finally, in order to ensure that the disks considered are dynami- 
cally stable, we only consider disks whose mass is lower than a 
tenth of the central star mass regardless of the value of ao- 

The observational determination of the disk lifetime orbiting 
stars of different mass could only be undertaken in recent years 
thanks to a growing number of observations. By correlating the 
fraction of stars showing an infrared excess with their mass for 
different stellar clusters, Kennedy and Kenyon (2009) were able 
to show that disks orbiting stars more massive than ~ 1.5Mq 
have shorter lifetimes. While a reliable and accurate dependency 
of disk lifetime is still difficult to extract from the current set of 
observations, the data seem to be consistent with disk lifetimes 

— 1/2 

decreasing as M^jJ , for Mstai- > I.SMq, and disk lifetimes in- 
dependent of the central star mass for less massive stars. In the 
following, we adopt this dependency but consider also a case in 
which the disk lifetime is independent of the stellar mass for all 
the masses considered (see section [373]i. 



10' 



3. Results 

3. 1 . Minimum mass versus semi-major axis diagrams 

Fig.|2]to|4]show the predicted mass versus semi-major axis dis- 
tributions of synthetic planet populations computed using either 
our nominal model (ffo - 1 .2), a model with a disk mass inde- 
pendent of the primary mass {ao = 0), and a model with a disk 
mass strongly dependent upon it (au - 2). In each of these fig- 
ures, the columns from left to right correspond to central stars 
with masses to 0.5, 1, and 2.0 Mq respectively. To allow for a 
comparison with observations, the middle and lower panels in 
each figures display the sub-population whose induced Doppler 
semi-amplitude is larger than 10 m/s or 1 m/s, and orbital period 
less than 5 year^ 

3.2. Tlieoretical mass function 

As planets migrate inwards during their formation, they might 
cross the inner boundary of our computational disk. When this 
happens, we freeze both the mass and the position of the planet. 
In our model, the final characteristics of these planets are there- 
fore ill-defined as they might be accreted by the central star, or 
might survive provided a stopping mechanism is at work in the 
inner parts of the disk. These planets might also grow further in 
mass due to accretion of gas flowing from the outer parts of the 
disk toward the inner regions. We do not take these effects into 
account and simply fix the position of these planets at 0.1 AU 
and keep the mass they had at the time they crossed this limit. 

Fig.|5]gives the planetary initial mass function for three val- 
ues of the primary mass (0.5, 1.0 and 2.0 Mq) and different val- 
ues of ckd- Three lines are displayed in each panel, corresponding 
to the full theoretical population and the planets that can be de- 
tected with a 10 m/s and a Im/s accuracy spectrograph (see cap- 
tion). The planets can be separated in three different classes (cor- 
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' Later, these planets will abusively be named "detectable with a 10 
m/s or 1 m/s accuracy spectrograph". 
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Fig. 2. Upper panels: mass versus semi-major axis diagrams for 
the nominal model. The mass of the primary is equal to 0.5 Mq 
(left panel), 1.0 Mq (central panel), and 2.0 Mq (right panel). 
Middle panels: mass versus semi-major axis diagrams for the 
nominal model, as can be detected with a 10 m/s accurate spec- 
trograph, and with a period lower than 5 years. The masses of 
the primary are the same as in upper panels. For each planet, the 
inclination angle between the line of sight and the orbital plane 
is chosen at random, in order to compute the resulting Doppler 
shift. Lower panels: as middle panels, but for a 1 m/s accurate 
spectrograph, and a period lower than 5 years. Note that in all the 
diagrams, planets crossing the inner boundary of the calculated 
disk are assumed to stay at 0.1 AU. The number of synthetic 
planets is equal to 30000 for each primary mass. 



responding to the three peaks clearly seen for Mjtar = 0.5Mq). 
The first class of planets are low mass planets (a few Earth 
masses). They correspond to planets having accreted all the plan- 
etesimals present in their feeding zone and have subsequently 
never grown massive enough to start type II migration. Because 
in the nominal model type I migration is strongly reduced, they 
remain close to their initial location and their growth is being 
quenched by exhaustion of the reservoir of planetesimals to ac- 
crete. Their final mass is therefore roughly given by the isola- 
tion mass at the initial location of the embryo. The second class 
of objects are Neptune-mass planets. These planets are formed 
in more massive disks, which allow them to become massive 
enough to start type II migration, thus increasing further their 
effective feeding zone. However, their mass never reaches the 
critical mass, hence they never enter runaway gas accretion and 
remain of Neptune-like structure. Finally, the third case corre- 
spond to planets equal or more massive than Jupiter, which form 
in quite massive disks and start their formation at large distance. 
These planets are able to reach the critical mass and enter a run- 
away gas accretion phase. The formation tracks of these planets 
(not shown here) are qualitatively very similar to the ones de- 
scribed in paper I, to which the reader is referred to for more 
details. 
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Fig. 3. Same as Fig. [2] but for a disk mass independent of the 
primary mass (q-d = 0). 
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Fig. 4. Same as Fig. [2] but for a-o - 2 (disk mass strongly de- 
pending on the primary mass). 



Our results show that the form of the mass function is sim- 
ilar for every value of ffo and primary mass. However, some 
differences can be seen. The most numerous planets are in the 
Super-Earth regime (below lOMf^), and their fraction does not 
strongly depend on the mass of the star, or the assumed an 
value. However, for - 2 around 2.OM0 stars, Super-Earth 
planets are slightly more massive (see the thickness of the first 
peak, which increases). Indeed, these planets essentially accrete 
all the mass in their feeding zone, reaching the local isolation 



Fig. 5. Planetary Initial Mass Function for different values of the 
primary mass (given at the top of each column) and different val- 
ues of ao (indicated on the right). In each panel, 30000 forma- 
tion models have been considered. The y axis presents the frac- 
tion of resulting planets in a given mass bin. The solid black line 
gives the mass function of all planets, whereas the red dashed 
and blue dotted lines correspond to planets observable assuming 
respectively a 10 m/s and 1 m/s accuracy spectrograph. 



mass. The isolation mass depends on the local surface density of 
solids (which increases for - 2 around 2.0Mq stars), and on 
the size of the feeding zone (which decreases for 2.OM0 stars, 
due to the increased gravity of the star). Therefore, planets are 
more massive only when the effect of the surface density is large 
enough to counteract the effect of the central star mass, namely 
for -2 and 2.0Mq stars. 

For higher mass planets, one can notice that the histograms 
for Neptune mass planets and massive planets evolves in an op- 
posite way. Neptune mass planets are more numerous around 
low mass stars compared to high mass stars. In particular, it is 
interesting to note that very few Neptune mass planets are pre- 
dicted around 2.0 Mq stars. This results essentially from the 
evolution of the disk mass (the effect is more pronounced for 
(Td = 2, see the maximum mass of planets around 2.0Mq), and 
can be explained as follows: when disks are more massive, plan- 
ets are more likely to accrete enough mass and to reach the crit- 
ical mass and enter in the "massive planets" population. 

Note that, as the primary mass increases, the size of the plan- 
ets feeding zone decreases (all the other parameters being kept 
constant), and one could conclude that this should result, on av- 
erage, in lower mass planets. However, this effect is clearly off- 
set by the increase in disk mass and the fact that, since many 
of the processes time scales vary with the Keplerian frequency, 
massive embryos form earlier around massive stars. All these 
results are qualitatively in agreement with Ida & Lin (2005) and 
Laughlin et al. (2004). While the formation of giant planets or- 
biting M dwarfs is indeed less probable, it is important to stress 
that such planets can form around M dwarfs within the core ac- 
cretion mechanism. 
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Interestingly enough, for primary stars ranging from 0. 1 to 
2.0 Mq the mass function for planets more massive than 100 Mg 
can be approximated by the following simple analytical law; 



10' I 1 1 1 r 



M, 



planet 



A^planet.lMo ^ 



Mo 



where Mpianet.iMo is the mass function around Solar mass stars, 
and 7 = 0.9 provides good results for the nominal case (ao - 
1.2) although the match is somewhat worse for stars below 0.3 
Mq. For a£, - 0, y - 0.8 provides a reasonable match, whereas 
for (Td = 2, larger values of y are required (around 1.3), but the 
match is significantly worse in particular for high mass stars. 

3.3. Close-in planets 

Massive planets orbiting at close distances to their stars (hot 
jupiters) are the easiest ones to detect by radial velocity surveys. 
So far, they represent about 10 - 20% of all the detected planets 
orbiting stars less massive than ~ I.SMq. Interestingly enough, 
for more massive stars jupiter-like planets are only discovered 
at significantly larger distances (greater than ~ 0.5 AU; see e.g. 
Wright et al. 2009). Since these planets are the easiest ones to 
detect, this lack of hot jupiters orbiting massive stars is really 
striking and cannot be accounted for by some observational bias. 

It has been proposed in the literature that this lack of mas- 
sive planets orbiting massive stars at close distances could be 
due to the shorter lifetimes of disks orbiting massive stars (see 
e.g. Currie 2009). In order to test this hypothesis, we have re- 
calculated the theoretical planet population orbiting 2.0 Mq stars 
but assuming that the distribution of disk lifetime is identical to 
the one orbiting solar type stars. Fig. [6] shows the comparison 
between the two resulting populations. We confirm the fact al- 
ready mentioned in the literature (e.g. Johnson et al. 2007) that 
assuming identical disk lifetimes results in the formation of a 
large number of hot planets orbiting 2.0 Mq stars which would 
have been easily detected. Conversely, reducing the liftetime 
of disks orbiting 2.0 Mq stars by 30 % shorter (following the 

— 1/2 

T'disk Ms,J suggested by Kennedy & Kenyon 2009 for high 
mass stars) results in a planet population without hot planets and 
only very few orbiting below 0.5 AU. 

This effect can be explained in terms of a reduced overall mi- 
gration for reduced disk lifetimes. Indeed, planets in orbit about 
massive stars grow larger in mass and more rapidly (see previous 
section) while at the same time the disk is evolving more rapidly. 
These two effects combine to hasten the transition from type II 
disk-dominated migration to type II planet-dominated migration 
which is the migration regime characterized by a ratio of planet 
mass to local disk mass exceeding unity. Compared to the nor- 
mal disk-dominated type II migration determined solely by the 
viscosity of the disk, the migration in this regime is significantly 
slower and eventually even comes to a halt as the ratio continues 
to increase. Hence, the overall extent of the migration is signif- 
icantly reduced and this despite the fact that higher mass stars 
are surrounded by higher mass disks, which, for a given planet 
mass, lead to more efficient migration in type I and type II disk 
dominated regimes. 

Note finally that for 0-0 = 2 and Mstai = 2.0Mq, a large 
number of warm super- Jupiter planets are predicted by the model 
(see Fig. |4|, even when the disk lifetime is reduced compared to 
disks orbiting solar mass stars. These planet are not detected by 
present day RV surveys, and therefore, such a value of au is 
excluded in the framework of our model. 
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Fig. 6. Histogram of final semi-major axis for detectable syn- 
thetic planets orbiting 2.0 Mq stars, assuming a 10 m/s accuracy 
spectrograph. Dotted line: disk lifetime are similar as for disks 
orbinting solar type stars, solid line: the disk lifetime is 30% 
shorter In each case, 9000 formation models have been calcu- 
lated. The number of observable planets is 424 and 621 respec- 
tively for the solid and dotted lines. 



3.4. Detection rate and properties of observed planets 

When comparing the planet population that formed orbiting stars 
with masses 0.5 Mq, 1.0 Mq and 2.OM0 in Fig.|2][3]and|4] we 
note that the number of detectable gas giant planets increases 
with the stellar mass considered. In all cases, we computed the 
formation of planets around 30' 000 stars and found respectively 
the following number of observable planets (with K > 10 m/s 
and a period lower than 5 years) 500-600, 900 and 800-1600 
depending upon the value of ao- 

We have compared the fraction of detectable planets for the 
different cases considered here, with the detection rate of planets 
orbiting different type of stars. Considering only planets more 
massive than five times Jupiter, and located between 0.5 and 2.5 
AU from the central star, Lovis and Mayor (2007) concluded 
that the detection rate is 5 times higher around A stars (5 planets 
meeting the criterion), compared to FGK stars. In the case of M 
stars, no such planets have been detected. However, it should be 
noted that their sample included stars much more massive than 
2.0 Mq. More recently. Bowler et al. (2009) have shown that the 
detection rate of similar planets around intermediate-mass stars 
is between 1.8 and 7 times higher than around solar type stars. 

We have compared the detection rate in our model for the 
0.5, 1.0 Mq and 2.0 Mq cases by normalizing the number of 
detectable planets orbiting stars of a given mass to the num- 
ber orbiting solar mass stars. In order to decide which planet 
is detectable, we used the same criterion as in Lovis and Mayor 
(2007). The ratio is found to range from 1 to 0.3 (for 0.5 Mq 
stars) and from 1 to 3 (for 2.0 Mq stars) for ckd ranging from 
to 2. Our ratio is therefore compatible with the two aforemen- 
tioned observational results only if large values of an are as- 
sumed. However, as pointed our in Sect. |3.3[ such a large value 
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Table 2. Results of Kolmogorov Smirnov tests comparing the- 
oretical and observed curves presented in Fig. |7] for the distri- 
bution of planetary mass (Mpianet) and planetary semi-major axis 

(flplanet)- 
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KS 


KS 


2.0 


1.2 


67 % 


76% 


2.0 





10% 


49% 


2.0 


2 


36% 


0.2% 


1.0 




0.2% 


0.3 % 



of ffD leads to large numbers of close-in massive planets orbiting 
2.0 Mq stars which are not observed. 

This apparent contradiction begs for an explanation. In fact, 
its origin could be three fold: Either the predicted masses are sys- 
tematically too low, or the predicted semi-major axis too large, 
or the absolute number of planets predicted to form is too low. 
As we shall argue below, we believe the later explanation is cor- 
rect and is a consequence of our assumption of a single embryo 
per system. 

We begin by looking at the cumulative histograms of min- 
imum masses and semi-major axis of the observed and syn- 
thetic planets orbiting stars of 0.5 and 2.0 Mq (Fig.|7]and[8]l. In 
these diagrams, we only considered planets with orbital period 
shorter than three years and inducing a Doppler semi-amplitude 
larger than a fixed value (see caption^ The actually detected 
planets meeting these criteria were taken from the Extrasolar 
planet encyclopaedia (http://exoplanet.eu/) as of December 1, 
2009. Finally, we plot on the same figures the distribution of the 
minimum masses and semi-major axis of the synthetic planets 
corresponding to the nominal model (a^, - 1 .2) assuming a 1 .0 
Mq Star 

From these figures as well as from the statistical tests pre- 
sented in Table |2] and Table |3] we first note that the nominal 
model calculated for solar type stars fails at reproducing the 
characteristics of the detected planets orbiting either M or A 
stars even though it was successful at reproducing the proper- 
ties of planets orbiting solar type stars (see paper II). 

As far as the planet population orbiting 2.OM0 stars are con- 
cerned, the masses and semi-major axis distributions obtained 
with the nominal model (q-q = 1 -2) are in reasonable agreement 
with observations. The main differences being that the models 
predicts slightly smaller semi-major axis and a larger number of 
massive planets (see Fig.|7]and Table|2]). 

In the case of planets orbiting M stars, the minimum mass is 
well reproduced for all values of ao considered while the semi- 
major axis distribution is not as well reproduced in particular 
for the a^, - 2 case. The main difference is an underestimation 
of the semi-major axis of planets located at distances larger than 
about 0.5 AU. However, note that comparing the semi-major axis 
distribution of synthetic planet orbiting 1 .0 Mq stars with the dis- 
tribution of actually detected planets orbiting M stars yields re- 
sults (see Table [3]) statistically considerably less significant thus 



Moreover, we note that only planets less massive than 20Mj at 
two sigmas are included in the Extrasolar planet encyclopedia (see 
http://exoplanet.eu/README.htmll. Therefore, it is difficult to com- 
pare models and observations in the very massive planets domain. We 
therefore exclude planets more massive than 25 Mj (to account in an 
approximate way for the 2 sigmas possible uncertainty in planetary 
masses) from the models and observations. 















/'/ //, 

y J A 





100 1000 10* 1 2 3 

U [Earth ma«s] « [AU] 

Fig. 7. Cumulative histograms of minimum mass (left) and semi- 
major axis (right) of synthetic and observable planets orbiting 
high mass stars. In each case, we only take into account planets 
whose induced Doppler semi-amplitude is larger than 30 m/s, 
and with periods lower than 3 years. Moreover, we exclude plan- 
ets located at distances lower than 0.1 AU (the internal limit 
of our computational disk). For synthetic planets (heavy lines), 
we use the formation model around 2.0 Mq stars, whereas we 
have selected planets discovered around stars with mass between 
I.8M0 and 3.OM0 (step-like lines). The three heavy lines are for 
the formation model around a 2.OM0 star, for cd equal to 1.2 
(solid line), (dotted line) and 2 (dashed line). The thin solid 
line is for a primary mass of 1.0 Mq. The results of statistical 
comparison between the different curves is indicated in Table [2] 

Table 3. Results of Kolmogorov Smirnov tests comparing theo- 
retical and observed curves presented in Fig. [8] 
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0.5 


1.2 


92% 


83 % 
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74 % 


78 % 
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87% 


55 % 


1.0 




11 % 


1 % 



showing again the importance of the role played by the mass of 
the central star. 

The fact that the properties of detectable synthetic planets 
orbiting both 0.5 and 2.0 Mq stars are similar to the observed 
ones lead us to believe that the discrepancy in the detection rate 
mentioned above is due to the fact that the number of plane- 
tary embryos increases with the mass of the central star From 
a purely logical point of view this has to be expected for two 
reasons. First, due to the scaling of the mass of the disk with the 
mass of the primary, more mass is available at a given metallic- 
ity to form embryos. Second, since these embryos grow faster in 
orbit of more massive stars, they are also more likely to become 
giant planets that will eventually be detected. 

Finally, we point out that for the planet populations orbiting 
both 2.0 Mq and 0.5 Mq stars, the nominal value for au provides 
systematically the best fit to observations. Interestingly enough, 
we recall that we determined the q-q = 1 .2 value not by fitting 
characteristics of planet population but by comparing accretion 
rates obtained by our alpha-disk model with observational deter- 
minations of this rate. Hence, amongst the values tested, the the 
nominal value of au is the value that consistently reproduce best 
both disk and planet observations. 



8 



Y. Alibert et al.: Extrasolar planet population synthesis. Ill 







1 






O.B 






Ml 


- 


B 

.2 0.6 
1 






....Mix 

"p " ' 




3 0.4 










0.2 




M j 










a. 






10 


100 


1000 


10* 




Fig. 8. Cumulative histograms of minimum mass (left) and semi- 
major axis (right) of synthetic and observable planets around M 
stars. In each case, we only take into account planets whose in- 
duced Doppler semi-amplitude is larger than 10 m/s, and with 
periods lower than 3 years. Moreover, we exclude planets lo- 
cated at distances lower than 0.1 AU (the internal limit of our 
computational disk). For synthetic planets (heavy lines), we use 
the formation model around 0.5 Mq stars, whereas we have se- 
lected planets discovered around stars with mass between 0.3Mq 
and 0.8Mq (step-like lines). The three heavy lines are for the 
formation model around a 0.5Mq star, for a£, equal to 1.2 (solid 
line), (dotted line) and 2 (dashed line). The thin solid line is for 
a primary mass of 1.0 Mq. The results of statistical comparison 
between the different curves is indicated in Table |3] 



3.5. Metallicity effect 

In our formation model, the metallicity of the central star enters 
only through the value of the planetesimals-to-gas ratio /d/g- 
The intrinsic difficulties in relating [Fe/H] and /d/g are discussed 
in Paper I, and we adopt here the same prescription, namely 
[Fe/H]= log(/D/G//D/G,G), with /d/g.g = 0.04 regai'dless of the 
mass of the primary (see Paper I for the justification of the /d/g,g 
adopted value). 

Fig.|9]shows the cumulative metallicity distributions of stars 
without planets and stars with a planet more massive than a given 
mass. To show the effect of varying the mass of the central stars, 
these distributions have been computed for stars with masses 0.5 
Mq, 1.0 Mq, and 2.0 Mq as well as for the three different values 
of au = 2. For reference, the metallicity distribution for all stars 
in the sample, regardless if they are hosting planets or not, is also 
drawn on the same figure. 

The overall effect of metallicity on the formation of plan- 
ets can be measured by the shift between the different colored 
curves (stars with planets) and the reference curve (all stars). The 
shift between the colored curves is a measure of the differential 
metallicity effect between stars of different masses but all hav- 
ing planets. We note that in all cases this differential metallicity 
effect is much smaller than the overall effect. Hence, we con- 
clude from our models that the metallicity effect is only weakly 
depended upon stellar mass. 

From the figure, it is clear that a metallicity effect is present 
in all cases in the sense that stars hosting planets are, on aver- 
age, more metal rich (colored lines are shifted with respect to 
the right of the black line). Comparing the three rows of Fig.|9] 
we also see that the metallicity effect also depends on the mass 
of the disk itself. Indeed, in the case of ao = the effect is 
very similar for the three stellar masses considered. On the other 
hand, for qto - 1-2 and, even more importantly for = 2, the 
effect increases as the mass of the primary decreases. This can 
be easily understood by recalling that the formation of a planet 



is, to first order, dependent upon the total amount of mass in 
form of planetesimals. In our model, this mass depends upon the 
product of the gas surface density (or the gas disk mass) and of 
the dust-to-gas ratio, which itself is a growing function of the 
metallicity. 

When the disk mass scales with stellar mass {au -1.2 and 
ffo = 2), the disks orbiting 2.0 Mq stars are much more massive 
than those orbiting 0.5 Mq stars. Therefore, the minimum metal- 
licity required to form a massive planet is correspondingly lower 
for massive stars than lower mass stars. Conversely, for = 0, 
the metallicity effect is much less dependent upon stellar mass 
since for all masses, the total amount of solids available is only 
varying with metallicity. 

Observationally, the metallicity effect is not defined by the 
existence of a planet more massive than a given value, but by the 
presence of an observable planet, with and induced Doppler shift 
larger than a given value. Fig. 10 present curves similar to Fig.|9] 
except that we now consider planets with periods less than five 
years but inducing a Doppler shift higher than 30, 10 or 1 m/s. 
In this case, both the mass and the semi-major axis of the planets 
are important as they enter both in determining the amplitude of 
the Doppler shift. 

For the two first cases, - and au - 1-2, the effect 
of metallicty on the distribution is somewhat higher for 2.0 M© 
stars. This is a consequence of the reduced migration resulting 
from shorter disk lifetimes already pointed out in Sect. 3.3 Since 
planets orbiting massive stars are formed, on average, at larger 
distances from their parent star, they induce a given Doppler shift 
only if they are massive enough. As seen above, larger planets 
are formed provided the total amount of solids is large enough. 
This condition is easier to fulfill in the case of a strong scaling of 
the disk mass with the stellar mass. It is therefore no surprise that 
the minimum metallicity required to form a massive planet (a 
planet inducing a large stellar reflex motion) is somewhat lower 
for larger mass stars when ao = 2. 



3.6. Composition 

The temperature in our proto-planetary disk model is computed 
as a result of the heat dissipation by viscosity, and therefore de- 
pends upon the gas surface density. As a consequence, the loca- 
tion of the iceline depends on the initial disk mass, as well as 
on the mass of the primary (viscous dissipation scales with the 
Keplerian frequency), as indicated in Eq. [T] We have compared 
the resulting composition of synthetic planets, for different val- 
ues of the au parameter, and for the three values of the mass 
of the central star considered so far The results are presented in 
Fig. 
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In order to differentiate between the bulk composition of dif- 
ferent planets, we have used different types of colors and sym- 
bols. Blue symbols represent planets that have accreted some 
icy planetesimals while red symbols are for planet having ac- 
creted only rocky planetesimals. The open symbols are for plan- 
ets where the mass of accreted gas is larger than the mass of ac- 
creted solids (envelope dominated) while the solid symbols are 
for planets that have accreted more solids than gas (core dom- 
inated). We note that this bulk composition corresponds to the 
time at which we end the calculation namely when the gaseous 
disk disappears. Later events that could modify this composition 
are not taken into account. Such events could be, for example, 
evaporation that could reduce the planetary gas mass envelope 
(e.g. Baraffe et al. 2004), or later accretion of icy bodies (see 
e.g. Morbidelli et al., 2000 in the case of the Earth). 
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Fig. 9. Cumulative histograms of metallicities of all stars and of stars harbouring a planet more massive than a given threshold (see 
top of the figure). The solid black line is the cumulative histogram of all stars considered in this work, which does not depend on 
the mass of the central star In each panel, the blue dotted curve is the cumulative histogram of 0.5 stars with a massive planet, 
the short dashed red line is for 1.0 stars, and the long dashed green curve is for 2.0 M© stars. 



From the figure, it is evident that the composition of our syn- 
thetic planets depends strongly on the mass of the primary star 
and, to a lower extent, to the value of a-o and hence the disk mass 
(the correlation between the disk properties and the planetary 
composition is studied in more details in Mordasini et al., 2010). 
Planets located closer than a few AU from their central star are 
found to be, on average, dryer when they orbit more massive 
stars. This is a consequence of the position of the iceline which 
varies in these different environments. Indeed, the iceline is lo- 
cated 1 to 2 AU further out when increasing the stellar mass from 
0.5 Mq to 2.0 Mq. This effect is further ampUfied for ao > 0, 



when the mean disk mass increases as a function of the primary 
mass. 



4. Discussion and conclusions 

We have extended the population synthesis calculations pre- 
sented in Mordasini et al. (2009a,b) to the formation of planets 
orbiting stars of different masses. Since the focus of this work 
was to pinpoint the effect of changing the mass of the central 
star and the corresponding changes in disk properties, we used 
the same model than presented in these two papers. In addition. 



10 



Y. Alibert et al.: Extrasolar planet population synthesis. Ill 

K=30 m/s K=10 m/s K=l m/s 



1 

0.8 



o 0.6 

as 



0.4 
0.2 
® 
0.8 



o 0.6 

Co 



0.4 
0.2 
(D 
0.8 



o 0.6 

Do 
In 



0.4 - 



0.2 




_ i 1 1 1 i 1 1 i 1 1 1 i i 1 i 11^. 

r // ~: 

- /'/ - 
r // i 

: y.^^' : 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 11^; 

// ^ 

/ // 

// ^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

/ i ^ 

1 1 1 r 'T'l^ 1 1 1 1 1 1 1 1 1 




1 1 1 1" 1 1 1 1 1 1 1 1 1 1 1 L-|-^r 

/ - 
/ --^ ^ 
~- 


/ '/ - 

// ^ 


1 fi 11 1 1 1 1 1 ' i^t^T: 
- //^' - 

//'^ ^ 
-^y^..t-^^\ , , , 1 , , , 1 ,- 


i^^: 

// ^ 

y.^¥---f\ , 1 , , , 1 , , , 1 ,- 


■Ti 1 rTi 1 1 1 1 1 1 1 1 1 11^ 
/f^ - 

-<r«Pl 1 1 1 1 1 1 1 1 1 1 1 1 



- o 



Q 
(0 

a 



(0 



C\2 



Q 
(0 

a 
"3 



II 

Q 

CO 

A 

a 



-0.4 -0.2 0.2 0.4 -0.4-0.2 0.2 0.4 -0.4 -0.2 0.2 0.4 

[Fe/H] 



Fig. 10. Cumulative histograms of metallicities of all stars and of stars with an observable planets. The solid black line is the 
cumulative histogram of all stars considered in this work, which does not depend on the mass of the central star. In each panel, the 
blue dotted curve is the cumulative histogram of 0.5 M© stars with an observed planet, the short dashed red line is for 1 .0 Mq stars, 
and the long dashed green curve is for 2.0 stars. A planet is observable if it produces a RV amplitude larger than a given K, 
and has a period lower than 5 years. The different assumed values of K and a-o are indicate on the top and right of the figure. Only 
planets more massive than the Earth are considered in the figure. 



we have explored the consequences on the planet population of 
difference scaling between the mass of the primary and the mass 
of the proto-planetary disk as well as between different relations 
between the lifetime of the disk and the central stellar mass. 

Our calculations show that the main effect of varying the cen- 
tral body's mass occurs through the variations in the disk prop- 
erties. The change in the later determine the changes in the syn- 
thetic planet characteristics. Assuming disk lifetimes indepen- 
dent of stellar mass results in the presence of many hot planets 



orbiting 2.0 stars, which is not consistent with observations. 
In addition, models assuming a very strong scaling of the disk 
mass with the primary mass (qto = 2) predict a lot of warm very 
massive planets orbiting 2.0 stars (see Fig. |4j right panels). 
This again seems to be excluded from present day planet obser- 
vations. As a consequence, we exclude from our model frame- 
work such high values of cd. On the other hand, if such high 
values were obtained from future disk observations, it would in- 
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Fig. 11. Composition of planets for different values of the primary mass and of otd- Red symbols are planets that have not accreted icy 
planetesimals, whereas blue symbols are for planet which have accreted a fraction of their solid mass in the form of icy planetesimals. 
Filled symbols are for planets whose envelope mass if lower than heavy elements mass, whereas open symbols are for planets whose 
envelope mass is larger 



dicate a that our models still lack some other physical mecha- 
nism. 

The dependance of the detection rate as a function of the 
primary mass is lower in our model as compared to observa- 
tionally inferred (Lovis & Mayor, 2007, Bowler et al. 2009). 
However, the predicted distributions of the planetary masses and 
semi-major axis around 0.5Mq and 2.0Mo stars are statistically 
comparable with the one observed around M dwarfs and A stars. 
We therefore conclude that the variation in the detection rate as 
a function of the central star mass are likely due to the fact that 
the number of planetary embryos, some of which leading to ob- 



servable planets, is a growing function of the primary mass. This 
in turn is likely to result from the larger amount of planetesimals 
available to form planetary embryos around massive stars. 

The first important result demonstrated in this paper is there- 
fore that the planetary population depends on the mass of the 
central star: indeed it is not possible to statistically reproduce 
properties of extrasolar planets discovered around low and high 
mass stars, using planet formation models developed for Solar 
type stars, as presented in Mordasini et al. (2009a,b). The sec- 
ond important result is that this dependance of planetary prop- 
erties on stellar mass is the result of changes in protoplanetary 
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disk properties for different types of stars. This shows the impor- 
tance of both future observational programs aiming at studying 
disk structure and properties, and the necessity to improve the 
description of protoplanetary disks in planet formation models. 
In particular, future papers will study the influence of irradiation 
by the central star (see Fouchet et al. 2010), and the presence of 
dead zones in these disks. 

Many future projects aim at detecting planets using tech- 
niques different from radial velocity surveys. One can cite tran- 
sit surveys (e.g. CoRoT and KEPLER), astrometry surveys 
(GAIA), microlensing, etc... A common point of all these sur- 
veys is that they target, or will target, stars of very different 
masses. This is for example the case for CoRoT and KEPLER 
which observe all stars in a given magnitude range and a given 
field. Comparing population synthesis models with data ob- 
served by these surveys is therefore intrinsically more compli- 
cated, since one needs to include also properties of the star-disk 
systems which are not yet well known. The results of these cal- 



culations are presented in Fig. 12 and can be used to compare 



theoretical results with the results of an observational survey. 
Such comparisons are presented elsewhere in the case of CoRoT 
(see Alibert et al. 2009), and for microlensing surveys (Alibert et 
al. 2010). As new observations, using different techniques, will 
become available, they will allow to put strong constraints on 
planet formation models. This will be the subject of future work. 
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Fig. 12. Mass versus semi-major axis for stars between 0.1 Mq and 2.0 Mq. The ao parameter is equal to 1.2, and disk lifetime are 
reduced for stars more massive than 1 .5 Mq. 30000 stars are considered in each panel. First line, from left to right, masses between 
0. 1 Mq and 0.5 Mq, second line, from left to right, masses between 0.6 M© and 1 .0 Mq, third line, from left to right, masses between 
1.1 Mq and 1.5 Mq fourth hne, from left to right, masses between 1.6 Mq and 2.0 Mq. 



